function powersA=computeAPowers(Hp, A)
oA = length(A);
powersA = zeros(Hp * oA, oA);
powerA = A;

powersA(1:oA, 1:end) = powerA;

for i=2:Hp
    [startPointer endPointer] = getPointers(i, oA);
    
    powerA = powerA * A;
    
    powersA(startPointer:endPointer, 1:end) = powerA;
end
end